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0 TURBULENCE 

o REACTION (treatment, kinetic schemes, emissions) 
o TURBULENCE/CHEMISTRY INTERACTIONS 
o ATOMIZATION 
o SPRAY EVAPORATION 

SIMULATION ISSUES: 
o NUMERICS (accuracy, convergence) 
o GEOMETRY (body-fitted grids, unstructured grids) 
o COMPUTATIONAL RESOURCES (Time, Storage) 


99 


PRECEDING 


PAG£ BLANK not FILMED 



JOINT VELOCITY-SCALAR PDF METHOD 


SIGNIFICANT MILESTONES AND RECENT PROGRESS 

o 2-D and 3-D time dependent flows (with finite-volume method) 

(Anand et al. 1987, Haworth & El Tahry 1989) 

o Stochastic dissipation model development and validation 
(Pope & Chen 1990, Pope 1991 , Anand et al. 1993) 

o 2-D Elliptic flows (mean pressure algorithm), swirling flows 
(Anand et. 1989, 1993) 

o Spray treatment 
(Anand 1 990) 

o Manifold methods for reaction kinetics 

(Maas & Pope 1992, 1994; Norris & Pope 1994; Norris & Hsu 1994) 


o Solve Poisson equation for mean pressure: 


3 2 <p> 

3x.dx. 
j J 


3 2 

3x.3x. 

i j 


<pU i U j > 


o Satisfy continuity by solving for velocity correction potential, velocity correction 


3 2 <j> 


= -ir <P u i> 


AU. = — 


36 

ox. 


5x.5x. 3 x l ' * * <p> 

o Solution algorithm is consistent with B-spline representation of mean fields 


o Same descretized form: A . s> = b 

o A is a banded matrix, constant 
and same for both <p> and <j> 

o LU decomposition only once 

o Special band solver economizes 
storage and computational effort 

o Judicious implementation of the 
algorithm results in significant 
economy in computer resource 
requirement 
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TURBULENT COMBUSTION MODELING ISSUES 


(FOR GAS TURBINE COMBUSTORS) 


o Most promising method for turbulent reacting flows 

ATTRIBUTES OF DIFFERENT PDF METHODS 


Method 

Joint PDF of <|> 


Joint PDF of 
U and ^ 


Joint PDF of 
U, |t, and (o 


Attributes 

Reaction treated 
exactly 


Reaction exact, 
Convection (mean and 
turbulent) exact, 
Variable-density effects 
exact 


Limitations/shortcominas 

Assumes gradient-diffusion, 

Does not give velocityfield 
(requires e.g, k-e) 
Turbulence/chemistry interactions 
not fully simulated 

Needs e equation 
(or equivalent) 


... In addition 

Provides complete closure, 
Treats turbulent streams of 
different scales, 

Can account for effects of 
large scale structures 


PDF CALCULATIONS FOR A RECIRCULATING FLOW 

(Anand et al. 1 989) 


'^r* f 


REATTACHHENT LEHCTH : 

EXPERIMENTAL (x„ t ) - 6.75 H 

CALCULATED - 6.65 H 


Backward-facing step 
Pronchick and Kline (1983) 


x * • 




0 3 a : a 3 


Storage: 1 .3 M words 
CPU Time: 6 min. Cray YMP 


<U>/U ref 
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STOCHASTIC DISSIPATION MODEL 


o Provides complete closure of the PDF equation (joint velocity-frequency-scalar) 

o More realistic than a mean dissipation model. Dissipation (rather, turbulent 
frequency) is also a random variable and included in the joint PDF. 

o Treats multiple scales in the flow 

o Accounts for internal intermittency 

o Accounts for effects of large scale structures, and influence of origin and history 
of the fluid particles 


dOO’ = -CO* <C0> (Su + C x Cl) dt + <(0> 2 h dt + CO* (2C X <00> cr 2 ) 1/2 dW 


dUj =-ii^dt + D i dt+(C 0 kco*) 1/2 dWj 
P d*i 


SWIRLING FLOWS 


o No theoretical limitations 

o Additional production terms due to non-zero mean swirl velocity 
o Additional terms in calculating the mean pressure (or mean pressure gradients) 

- Boundary layer flows: 

> radial pressure gradient 

> axial pressure gradient also included 

- Elliptic flows 

> additional terms in the Poisson equation for pressure 

o Validation of the stochastic dissipation model and first calculation of swirling 
flows with the joint PDF method (Anand et al. 1993) 
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JOINT PDF CALCULATIONS FOR SWIRLING FLOWS 


COMPARISON WITH REYNOLDS-STRESS MODEL RESULTS AND 
ASSESSMENT OF GRADIENT DIFFUSION MODELING 


x/D = 1 .06 2.65 5.29 



O.H »— J 1 f— 

-.40 -.20 0.0 .20 .40 


1 000<u2v>/<U>3qc 



PDF CALCULATIONS 

RS MODEL WITH RSM RESULTS 

RS MODEL WITH PDF RESULTS 
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SPRAY CALCULATIONS 

(Anand 1990) 


o Advanced spray models (stochastic 
lagrangian, Monte Carlo) naturally 
compatible with the joint PDF method 

o Assumptions about turbulent kinetic 
energy partition avoided 

o Effects of gas phase turbulence 
structure (velocity cross-correlation) 
included 


■ ''0*0.62 1.04 2.07 3. H 6.22 12.4 



Computed profiles of normalized 
turbulent kinetic energy of air compared 
against data. 


105 micron glass beads. NASA HOST C data 


w/0- a. 62 I . 


2 .07 3 . i t 



<Ud>A^o.c 


«/!>• 0.62 1.04 2.07 3,il 6.22 12,4 



10 U'(j/U d O.C 


REDUCED KINETICS / MANIFOLD METHODS 


o Low dimensional manifold methods (ILDM, TGLDM) 

- Given detailed kinetes, they provide low-dimensional description 
(e.g., 1-D, 2-D, 3-D) in multidimensional composition/scalar space 

- Use dynmical systems theory to determine the low. dim. manifold 

- Avoid ad hoc assumptions, e.g, partial equilibrium of some of the reactions 
Implications for ignition and lean blow-off 

- Not fuel specific like conventional reduced kinetic schemes 


0.006 r~ 



Perfectly Stirred Reactor (Pope & Maas 1993) 



'CH 1 "H 

5.0 10° 


J 2.0 1C' 3 


1.0 io ,J 


Laminar Premixed Flame (Maas & Pope 1994) 
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PARALLEL PROCESSING 


o Objective: Turnaround time of Iday or less for 3-D combustor calculations 
o Particle partitioning, domain decomposition (multigrid, multi-block) 

o Preliminary results for 2-D flow with particle partitioning (Pope 1994) 

- 16 nodes, 128 MB each, IBM SP1 

- 12.8 million particles (800,000 per processor) 

- 50 time steps 

- 44 minutes/processor (45 minutes clock time) 

Extrapolation to 3-D combustor calculations 

- 6.5 hours clock time with 32 processor SP1 


JOINT PDF FOCUS AREAS 


o 3-D Flows, Improved solution algorithms 
o Parallel processing 

o Reduced kinetics / Low Dimensional Manifolds 
o Evaporating / reacting sprays 
o Emphasis on emissions and performance predictions 
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